Ideal glass transition in a simple 2D lattice model 
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We present a simple lattice model showing a glassy behavior. R matrix analysis predicts critical 
termination of the super-cooled fluid branch at density pg = 0.1717. This prediction is confirmed by 
dynamical numerical simulations, showing power-law divergences of relaxation time ri/2, as well as 
the 4-susceptibility X4 peak's location and height exactly at the predicted density. The power-law 
divergence of X4 continues up to X4 as high as 10'*. Finite-size scaling study reveals divergence of 
correlation length accompanying the transition. 



Understanding the transition of supercooled liquids 
into a glass is considered by many to be one of the out- 
standing challenges of condensed matter physics. Many 
liquids, when cooled fast enough to avoid crystallization, 
appear to freeze into solid-like structures devoid of crys- 
talline order [J |3) The time scales for structural 
relaxation in such metastable super-cooled regimes in- 
crease dramatically as temperature is lowered. For strong 
glass formers the relaxation times grow exponentially 
T = Toexp(A/T). Fragile glass formers exhibit relax- 
ation times that increase more rapidly than Arrhenius 
and are often fitted by Vogel-Tammann-Fulcher (VTF) 
functional form t = roexp[yl/(T — Tq)], with a charac- 
teristic temperature Tq !4| . 

The glass transition temperature is experimentally de- 
fined as the temperature at which dynamic relaxation 
times exceed those accessible in typical experiments, e.g. 
when viscosity hits 10^^ Poise. During the past century, 
much interest has been focused on understanding the na- 
ture of this transition. Clearly, the glass transition tem- 
perature, defined by some viscosity cutoff value, is just 
an arbitrary reference point along the gradual increase 
of relaxation times with decreasing temperature. The 
question whether there is some deeper, profound, physi- 
cal meaning to the glass transition is still hotly debated 
[l[ . Is the fast increase of relaxation times merely a sharp 
cross-over in the dynamics, or could it be manifestation 
of a true thermodynamic transition? [Obviously, when 
considering a glass transition in a system exhibiting a 
solid phase, such as real glass formers, the notion of a 
thermodynamic glass transition must be interpreted in 
the sense of a restricted part of phase space. For sim- 
plicity, we ignore this distinction in the following.] Many 
theoretical studies have been applied to support either 
one of the competing views. For example, a popular 
microscopic approach is the mode-coupling theory (see 

Q for reviews). It predicts a dynamic glass transi- 
tion, characterized by ergodicity breaking, while ther- 
modynamic (equilibrium) quantities such as the isother- 
mal compressibility do not become singular. In contrast, 
the replica approach [tII predicts a structural glass tran- 
sition with pure thermodynamic origin, characterized by 
a vanishing configurational entropy. Other phenomeno- 
logical theories, such as the random first order transi- 



tion (RFOT)!, and the potential energy landscape 
(PEL)|l3| to name only two, also predict a thermody- 
namic phase transition. 

Recently, the R matrix ^11', IS] approach for analysis 
of the Mayer cluster inte gral s expansion has been applied 
to the hard-spheres fluid IJ, llj| . It provides the density 



as a function of the activity z {z = e^'^, where ^ is the 
chemical potential and [3 is the inverse temperature) and 
predicts a critical termination of the super-cooled fluid, 
with a power-law divergence of the isothermal compress- 
ibility. The packing-fraction at which this divergence is 
predicted to happen is 0.556(5), surprisingly close to the 
experimentally reported glass transition packing-fraction 
0.56(1). This result, therefore, strongly supports the 
existence of a thermodynamic glass transition for hard 
spheres underlying the (experimentally and numerically 
observed) dynamical arrest. It is desirable to have nu- 
merical measurements of the super-cooled hard-spheres 
equation of state near the transition in order to test the 
validity of the R matrix approach. However, these sim- 
ulations are extremely challenging. Accordingly, contra- 
dicting results have been reported regarding the existence 
of singularities in thermodynamic quantities for this sys- 
tem 



The limits of numerical methods often hamper the 
study of glass transition. Excluding the non-physical ki- 
netically constrained models, most models studied are 
either complex (binary mixtures) or hard to simulate 
(hard spheres). They are therefore limited in system size 
and simulation times. For example, a recent study of 
Lennard-Jones binary mixture Q reports that enlarg- 
ing the system to include 27000 particles improves the 
quality of the extrapolation of /c-dependent quantities to 
zero wave vector. Moreover, simulations are generally 
limited to time scales roughly ten orders of magnitude 
shorter than those near the laboratory glass transition 
temperature Tg and therefore to the initial stages of the 
glass formation process 18| . These numerical limitations 
might be lifted by introducing a simpler model system 
that still captures the essence of glassy behavior. Keep- 
ing that in mind we set to explore the glass transition is 
the iV3 lattice model. 

The iV3 model is a simple 2D model on a square lat- 
tice. Particles interact only through hard-core exclusion 
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up to the S""*^ nearest neighbor. The model is known to 
undergo a first order solidification transition 1^ 2^, 2ll 
where density jumps from p/ — 0.161 to Ps — 0.191 [21[ 
(the closest packing density is 0.2). Like the hard spheres 
case, R matrix analysis predicts a critical termination of 
the super-cooled fluid where the isothermal compressibil- 
ity power-law diverges. The critical density is found to 
be pt ~ 0.1717. In concordance with hard spheres results 
[ij] , we hypothesize that this point is indeed the thermo- 
dynamic glass transition for this system. We then study 
the dynamics of the model by extensive MC simulations 
and find that the dynamical quantities diverge exactly at 
the density predicted. We therefore conclude that the dy- 
namical arrest in the iV3 model results from a singularity 
of the free energy, as predicted by the R matrix. These 
results support the view of a thermodynamic (a.k.a ideal) 
glass transition in this system. Furthermore, we propose 
the A''3 system as a simple and convenient model-system 
for future studies of glassiness. 

In order to construct the R matrix for the iV3 model, 
we extended the number of known Mayer cluster inte- 
grals to 23, using the transfer matrix (TM) method. 
We have employed a diagonal-to-diagonal, symmetry re- 
duced, TM, with strip width as large as M = 24 (3874112 
symmetry-reduced classes). The cluster integrals pro- 
vide the exact 11x11 leading R submatrix presented in 
table H] (for details on R matrix construction, see [l^). 
The matrix elements quickly converge to a well-defined 
asymptotic form which we use to extrapolate additional 
matrix elements and obtain the equation of state (figure 
[1]). Remarkably, the results, based only on low-density 
expansion, are in an excellent agreement with both MC 
data and exact TM calculations. The physical singularity 
is found at zt ~ 66.67, well above the first order transi- 
tion [zc — 39.496), with a critical density pt = 0.1717. 
Furthermore, the R matrix provides an exact formula for 
the critical exponent a' associated with the termination 
point of the fluid [13, 14 1 : near this singularity, the den- 
sity is given by 



Pt - p{z) ~ {zt - zf , 



(1) 



and the critical exponent is found to be a' = 0.39(2). 

We hypothesize that this thermodynamic criticality 
underlies an ideal glass transition for the N3 model, and 
set out to study the model dynamically looking for sig- 
natures of this glass transition. We conducted canonical 
(constant density) MC studies of the model in the fol- 
lowing way: The starting configuration was generated 
under extreme cooling conditions (or, equivalently, infi- 
nite chemical potential). Particles were allowed to diffuse 
when no insertion was available. This process is known 
to terminate at the random closest packing (RCP) state 
with density 85% of closest packing density [i^]. Here, 
we stop the cooling at the desired density (below RCP), 
and let the system relax diffusionally. Given enough time, 
the global equilibrium phase-separation state is reached. 
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FIG. 1: (Color Online) N3 equation of state: R matrix pre- 
diction, based on the first 23 Mayer cluster integrals (dashed 
line), Monte-Carlo calculation on a 1000 x 1000 lattice (solid 
line), and exact transfer matrix calculation for a semi-infinite, 
25 sites wide, strip (symbols). The latter two methods provide 
equilibrium results, while the R matrix extrapolates to the 
super-cooled fluid branch. The agreement of the R matrix re- 
sults with the numerical methods is excellent throughout the 
fluid regime. Inset shows the diagonal (-B,i) and ofi'-diagonal 
{An) R matrix elements, together with the fitted asymptotic 
form. 



On shorter time scale the system relaxes to a disordered 
phase. We first measure the density-density correlation 



C(t) = ^ (1 / <n(r,0)n(r,i)>dr-p| (2) 



along the relaxation process (n(r, i) = 1 if a particle 
exists at site r at time t and otherwise). Figure [2] 
shows the typical glassy dynamics picture: a plateau 
(/3 regime) followed by a stretched-exponential decay (a 
regime). Due to the discrete nature of the diffusion pro- 
cess in this model, the /3 relaxation stage is very short 
(of order one simulation time unit) and is not presented. 
The relaxation time ti/2, defined as the time at which 
Cii) — 1/2, power-law diverges as the density approaches 
Pg = 0.1717: 



Tl/2(P) ~ (pg -P) ^ 



(3) 



with p, — 0.83, (figure [2l in set). In addition, we measure 
the 4-susceptibility ya [23l. [ii ] 



X4(i) - A^(< C(tf >-< C{t) >2). 



(4) 



Again, a typical glassy behavior is observed (figure [3|) - 
X4 peaks at the a phase, and the peak grows in height and 
shifts to higher times as density increases. Peak heights 
iXmax) and locations (T4) also power-law diverge as pg 
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TABLE I; Mayer cluster coefficients nbn and R matrix diagonal {B„) and off diagonal {An) elements for the N3 model 




FIG. 2: (Color Online) Density-density auto-correlation ((2|, 
showing the typical glassy relaxation picture, on a logarith- 
mic time scale. The inset shows the relaxation time ri/2 vs. 
density (symbols), which is well-fitted by a power-law ([3]) di- 
verging at pg. 

is approached (figure [3l inset). A transition to activated 
dynamics occuring close to the glass transition could be 
manifested by the onset of a slower, logarithmic, growth 
of the Xi peak [25| . We do not observe any such transi- 
tion for X4 values up to 10'*. 

The above MC data confirm the R matrix prediction 
to an excellent agreement. Given that this prediction is 
based solely on low-density series expansion, it is remark- 



FIG. 3: (Color Online) X4 a-s a function of time showing the 
familiar peaks. Peaks' heights, Xmax, and locations, Ti, power- 
law diverge as density approaches pg (inset) with critical ex- 
ponents 2.85 and 0.93 respectively. 

able that it captured quantitatively the behavior at the 
deep super-cooled regime. This attests for the validity 
of the R matrix approach and its prediction of a ther- 
modynamic criticality in the equation of state of the N3 
super cooled fluid, and provides a strong evidence that 
the glass transition in this model is indeed a thermody- 
namic, ideal, one. 

The growing Xi peak is indicative of growing coop- 
erative correlations in the relaxation process [26l |. It 
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FIG. 4: (Color Online) Finite size analysis: Xi peak height as 
a function of lattice linear size, for various densities. As the 
density approaches pg, larger lattices are needed for converged 
results, attesting for a diverging length scale. For lattices 
smaller than the correlation length, Xmax is expected to grow 
like system's size. A straight line with a slope of 2 is presented, 
to guide the eye. 



measures the volume upon which difFusional moves are 
correlated 27| . In concordance, growing correlations 
lengths are seen also by the emergence of finite size ef- 
fects in the density-density correlations as shown in fig- 
ure m These finite-size effects, recently highlighted by 
Karmakar et al. 28] underscore the importance of using 
large systems for MC studies of glassiness, which is most 
difficult for popular models currently used. 

Unlike the hard spheres case, the termination density 
Pg predicted by the R matrix and the dynamical arrest 
occur in close proximity to the random closest packing 
density prcp [22] ■ Therefore, the region beyond the tran- 
sition is inaccessible in this model. It is important to 
note that the co-occurrence of the two phenomena is not 
a universal trait of the i?-matrix analysis. For example, 
in the hard-spheres model the R matrix prediction for 
the ideal glass density is pt = 0.556(5), much lower than 
Prcp — 0.64. Hopefully, future work will find a model 
that is as simple as the N3 but also allows access to den- 
sities beyond pt ■ This could be achieved by studying the 
soft-core N3 model, or other hard-core lattice models. 

We stress that the simplicity of the N3 model is impor- 
tant not only in order to allow for analytical treatment, 
but to facilitate numerical studies of large systems, much 
larger than those typically used in glass studies. This is 
especially important when one approaches the glass tran- 
sition, where long-range cooperative relaxation processes 
emerge, manifested by significant finite-size dependence. 
For example, at density p = 0.1715, even a 1000 x 1000 
lattice (171500 particles; linear size ~ 447 particle diam- 
eters, much larger than typical 3D studies) is not large 



enough to converge to bulk values as seen in figurelH The 
need for a simple model then is not a matter of comfort 
but a real necessity. We therefore propose that the N3 
model, or similar models, could serve in future studies of 
glass formers being simple to handle, yet capturing the 
essence of glassiness. 

In conclusion, we have applied the R matrix approach 
to the N3 model and found that its super-cooled equa- 
tion of state becomes singular at density pg = 0.1717, 
where the isothermal compressibility power-law diverges. 
MC simulations confirm that the model shows the char- 
acteristics of a fragile glass former undergoing a glass 
transition at the predicted pg . It thus follows that in this 
model the phenomenological glass transition, observed as 
a fragile-glass dynamical arrest at pg, is accompanied by 
a thermodynamic criticality. 
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